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Abstract 

One of the main drawbacks of standard cosmology, known as the horizon 
problem, was until now thought to be only solvable in an inflationary sce- 
nario. A delayed Big-Bang in an inhomogeneous universe is shown to solve 
this problem while leaving unimpaired the main successful features of the 
standard model. 
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1 Introduction 

Standard cosmology is known to rest on three observational pillars: the expansion of 
the universe following Hubble law, the nearly isotropic black body cosmic microwave 
background radiation (CMBR) and the abundances of light elements produced dur- 
ing nucleosynthesis. 

Besides these successfull predictions, it leaves ununderstood other peculiar fea- 
tures of the observed universe. 

In the present letter, a large class of initial singularity surfaces, the study of 
which has been initiated in a previous work [|]], will be used to address one of the 
drawbacks of standard cosmology: the horizon problem. 

The problem is the following: in hot Big Bang (BB) universes, the comoving 
region over which the CMBR is observed to be homogeneous to better than one 
part in 10 5 is much larger than the comoving future light cone from the BB to the 
last scattering surface. The latter provides the maximal distance over which causal 
processes could have propagated since a given point on the BB surface. Hence, the 
observed quasi-isotropy of the CMBR remains unexplained. 

Solving this problem was one of the main purposes of the inflationary paradigm as 
it was first put forward by Brout, Englert and Gunzig || in 1979 and independently 
by Guth || in 1981. But, as inflationary scenarii inflatened, some self-produced 
undesirable features came into the way: for instance, reheating is not actually well 
understood and important details are still under study j4|. 

But despite these drawbacks, inflation has by now become a quasi- standard 
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paradigm as it was thought to be the only way to deal with the major horizon 
problem. 

Following a suggestion by Hu, Turner and Weinberg ||, Liddle @ has even 
proposed a proof that inflation is the only possible causal mechanism capable of 
generating density perturbations on scales well in excess of the Hubble radius, and 
hence the only way of solving the horizon problem. 

As it was stressed by these authors, this problem involves the homogeneity and 
isotropy of the Freedmann-Robertson- Walker (FRW) model, proceeding from the 
Cosmological Principle upon which rests standard cosmology. Liddle's entire argu- 
ment depends only on the properties of the FRW metric. 

The so-called Cosmological Principle is in fact not an a priori principle, but at 
most a simplifying working hypothesis: the universe being as it is, all astrophysics 
can do is to build models compatible with observation, should they contradict the 
Cosmological Principle. We come back to this point in the discussion. 

The purpose of this letter is to solve the horizon problem by means of a delayed 
BB singularity in an inhomogeneous model of universe, thus discarding this Princi- 
ple. For simplicity, we use a Tolman-Bondi model, since it allows a fully analytical 
exact reasoning. 

This model will be described in section 2. Calculations and arguments will be 
developed in section 3 and some examples given in section 4. Section 5 will be 
devoted to a brief discussion of the results and to the conclusion. 
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2 An inhomogeneous delayed Big Bang model 

In an expanding universe, going backward along the parameter called the cosmic 
time t means going to growing energy densities and temperatures. 

As one goes down the past, from our present matter dominated age defined 
by the constant temperature hypersurface T ~ 2.73°K, one reachs an epoch when 
the radiation energy density overcomes the matter one. This radiation dominated 
area lasts until Planck time, Tpi ~ 10 19 GeV, which marks the limit beyond which 
quantum gravitational effects are expected to confuse our understanding of the laws 
of physics. 

To deal with the horizon problem, one has to compute light cones. As will be 
further shown, a large class of models can be found for which the horizon problem is 
solved by means of light cones never leaving the matter dominated area. We will thus 
retain the Tolman-Bondi model for dust, an ideal non zero rest mass pressureless 
gas. 

2.1 The class of Tolman-Bondi models retained 

The Bondi line-element M, in comoving coordinates (r,9,tp) and proper time t, is: 



Solving Einstein's equation for this metric with the dust stress-energy tensor 
gives: 



ds 1 = -c 2 dt 2 + S 2 (r, t)dr 2 + R 2 (r, t){d0 2 + sin 2 Odip 2 ) 



(1) 



S 2 (r,t) 



R' 2 (r,t) 



(2) 



1 + 2E(r)/c 2 
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47rp(r, t) 




(4) 



(3) 



R'(r,t)R 2 (r,t) 



where a dot denotes differentiation with respect to t and a prime with respect 
to r. p(r,t) is the energy density of the matter. 

E(r) and M(r) are arbitrary functions of r. E(r) can be interpreted as the total 
energy per unit mass and M(r) as the mass within the sphere of comoving radial 
coordinate r. 

M(r) remaining constant with time, it is used to define a radial coordinate r: 
M(r) = M r 3 , where M is a constant. 

Equation (|3]) can be solved and gives a parametric expression for R(r, t) for 
E{r) 7^ and an analytic one for E(r) = 0. 

As there are evidences that the observed universe does not present appreciable 
spatial curvature, it can be reliably approximated by a flat E{r) = Tolman-Bondi 
model. 

With the above definition for the radial coordinate r, R(r,t) possesses thus an 
analytical expression, which we write: 



to(^) is another arbitrary function of r, representing the BB singularity surface 
for which R(r,t) = 0. One can always choose t (r) = at the center (r = 0) of the 
universe by an appropriate translation of the t = const, surfaces and describe our 
universe by the t > t (r) part of the (r, t) plane, increasing t corresponding to going 
from the past to the future. 
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Equation @ substituted into equation (|j) gives: 

1 

P[r,t) ~ 27rG[3t-3t Q (r)-2rt' (r)][t-t (r)] 

2.2 Shell-crossing 

The above expression for p leads to two undesirable consequences 

1) The energy density goes to infinity not only on the BB surface t = to(r), but 
also on the shell- crossing surface: 

t = t (r) + ~rt' (r) (7) 

2) This energy density presents negative values in the region of the universe 
located between the shell-crossing surface (|7]) and the BB singularity, corresponding 
to 3t — 3to(r) — 2rt' (r) < and t — t (r) > 0. One can wonder what does physically 
mean a negative energy density for dust. 

Shell-crossing is thus generally considered as a mischief of Tolman-Bondi models 
and physicits usually try to avoid it ||, e.g. by assuming t' Q (r) < for all r. 

But, as will be developed in next section, we need an increasing BB function 
£o(V) to solve the horizon problem. Let us hence briefly show how to circumvent 
these two difficulties while keeping t' (r) > 0. 

1) A way out the shell- crossing surface problem is to consider that, as the energy 
density increases while reaching its neighbourhood from higher values of t, radiation 
becomes the dominant component of the universe, pressure can no more be neglected 
and the Tolman-Bondi model does no longer hold. 
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(6) 



2) A negative value of p proceeds from a negative value of R' in equation (f|) . 
The physical definition of energy density is: 

SM being the element of mass in an element of volume SV. 
The element of 3-volume corresponding to the flat Tolman-Bondi metric, i.e. 
metric (0) with S 2 (r,t) = R' 2 (r,t), is: 

6V = R'R 2 sin 6drd6d V (9) 

which, when integrated over 6 and (p, becomes: 

5V = 4irR'R 2 dr (10) 

As the physical volume 5V is by convention always positive, equation ( |T0| ) pos- 
sesses a physical meaning only if it is written: 

5V = An\R'\R 2 dr (11) 

And thus, in equation (|J), one has to replace R' by \R'\, which gives in equation 

P(r ' ^ = 2nG\3t - 3t (r) - 2rf (r) \ [t - t (r)] ^ 
However, as, in the following, the light cones of interest never leave the region 
situated above the shell-crossing surface in the (r, t) plane, 3t — 3t (r) — 2rt' Q (r) 
remains positive and equation @ holds. 
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2.3 Definition of the temperature 

In the course of this letter, we shall be led to use surfaces of constant tempera- 
ture T. Since the universe is not homogeneous, there is, at a given t, no global 
thermodynamical equilibrium, and T is not readily defined. We assume that the 
characteristic scale of the p inhomogeneity is much larger than the characteristic 
length of the photon-baryon interaction and that there is always a local thermody- 
namical equilibrium. This enables us to define a local specific entropy S by: 

o (r) _ k B n^r,t)m b 

where m& is the baryon mass and k B the Boltzmann constant. 
We then define T by: 

n 7 = a„T 3 (14) 

where a -- >" /. 



n 7r 2 (7k:) 3 

The following expression for T can then be obtained from equations (|12|) to ([14|) : 

T ( rt \ = ( ?W V /3 ( 15 \ 

K,J \27iGk B m b a n \3t - 3t (r) - 2rf (r)\[t - t (r)] J 1 ; 

The equation of the T = const, surfaces located after the shell- crossing surface, 
i.e. with 3t — 3to(r) — 2rt' Q (r) > 0, is thus the positive solution of the second order 
in t equation derived from equation (Jl5|): 
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2.4 The "centered Earth" assumption 

In this first approach of a delayed BB solution, the Earth will be assumed situated 
sufficiently close to the "center" of the universe, so as to justify the approximation 
r p = 0, the subscript p refering to our actual location at the present time. We shall 
comment on this "center" of the universe in the final discussion. 

The value S p of the entropy function at (r p , t p ) is S p = kBi] p , r\ p being the present 
local photon to baryon density ratio, which is taken to be of order 10 8 . 

Remembering that t (r=0) has been chosen to be zero, one can add to the 
specifications of the t (r) function: 

rt[,| r= o = 

to get at (r = 0): 

t p ~ 3.10 17 s for T p = 2.73° K 

t is ~ 6.10 12 s for T ls = 4000°fT 

These values are of the same order of magnitude as in the standard hot BB 
model. The nucleosynthesis scenario would thus approximately be the standard one 
for the here described universe in the vinicity of r = 0, provided the characteristic 
lenght of the density inhomogeneities is much larger than the mean free path of the 
nucleons. This latter condition will be discussed in section 5. 

As the light elements abundances predicted by standard cosmology fit rather well 
the data observed in our neighbourhood, the choice of the "center" of the universe 
for the location of the Earth seems justified, as far as the above cited condition 
obtains. 
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This good agreement between the observed data and the abundances predicted by 
the standard model led Liddle |J to adopt the Cosmological Principle and thus assert 
that the FRW metric obtains for the whole universe. This was a key-assumption for 
his tentative proof discussed in our Introduction. 

But this assumption seems exceedingly narrowing as far as the present available 
data have been measured in our direct neighbourhood as compared to cosmological 
distances - the today most remote measured abundances are for deuterium at red- 
shifts z < 5 ||. The "centered Earth" assumption seems thus enough to complete 
the game. 

3 Solving the horizon problem 

Light travels from the last scattering surface to a present local observer on a light 
cone going from (r p = 0, t p ) to a 2-sphere (ri s ,ti s ) on the last scattering 3-sphere 
defined by T = 4000° K. 

To solve the horizon problem, it is sufficient to show that this 2-sphere can be 
contained inside the future light cone of any (r = 0, t > 0) point of space-time. 

One of the key-points of the reasoning here proposed is a shell-crossing surface 
situated above the BB surface and monotonously increasing with increasing r, which 
is always verified if: 

t' (r) > for all r 
5t' (r) + 2rtl(r) > for all r 

A to( r ) function increasing with r implies that the BB "occured" at later t 
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for larger r, hence the evocative "delayed Big-Bang" we choosed to qualify this 
singularity. 

The model being spherically symmetrical and the null geodesies being radial, the 
relevant light cones are obtained for 9 = const, and ip = const, in equation ([!]). 
Writing ds 2 = in this equation, one gets a differential equation for the null 



cones: 



# = ±- (17) 

dr c 

Substituting above the expression of R' obtained from equation (|5]), one finds: 



dt , 1 /9GM \ 1/3 3* - 3t n (r) -2rt 



'o(r) 



dr ± 3c{ 2 ) [t-to(r)] 1 / 3 
Comparing to equation ([7|), one immediately sees that the curves representing 
the light cones in the (r, t) plane possess an horizontal tangent on and only on the 
shell-crossing surface, where 3t — 3to(r) — 2rt' (r) goes to zero. 



The curve t(r) for the past light cone from (r p , t p ) verifies equation fll8|) with the 
minus sign. 

As far as one considers the part of this light cone located after the shell-crossing 
surface, and thus after the BB singularity, 3t — 3to(r) — 2rt' (r) and t — t (r) remain 
positive and is always negative. t(r) is a strictly decreasing function of r and the 
light cone will have to cross the strictly increasing shell-crossing surface at a finite 
point where the derivative of t(r) goes to zero. 

On its way to shell-crossing, the null geodesic will cross in turn each T = const, 
surface at a finite point. 

Let (ri s i,ti s i) be the coordinates of the crossing point on the last scattering 
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surface T = 4000° K. 

Now consider a backward null radial geodesic starting from any point above 
the shell-crossing surface and directed towards r = 0. Its equation is a solution of 
differential equation ([TJ) with the plus sign. 

Its derivative remains positive as long as it does not reach the shell-crossing 
surface. If it was to reach this surface before the "center" of the universe, its 
derivative would go directly from a positive value to zero, which would imply for 
the curve of the light cone an horizontal tangent in the (r, t) plane. 

Since we consider models for which the curve representing the shell-crossing sur- 
face is strictly increasing with r, it cannot be horizontally crossed from upper values 
of r and t by a strictly increasing curve. And thus one is led to an inconsistency. 

This implies that the backward light cone starting from any point above the 
shell-crossing surface reachs r = at t c without crossing this surface, and thus with 
t c > 0. 

This statement holds for every light cone issued from any point on the last 
scattering surface. 

There is thus an infinite number of points (r = 0, t c > 0) of which the future light 
cone contains the sphere on the last scattering surface seen today in the CMBR. 

Every point on this sphere can be causally connected and the horizon problem 
is solved. 
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4 Examples of appropriate Big Bang functions 

In previous sections, conditions have been imposed upon the BB function t (r). 
They can be summarized as follows: 

t o (r = 0) = 

*o(r) > for all r 

5t / (r) + 2rt" (r) > for all r 

rto| r=0 = 

It is easy to verify that the class of functions: 

t (r) = 6r n 6 > n > (19) 

fulfills these conditions. 

Another feature imposed upon the model to justify the dust approximation is 
that the light cones, for the Tolman-Bondi metric, never leave the matter dominated 
area. This prescription has been tested upon peculiar models of the above class with 
S(r) = const. = kBi] p - 

A number of light cones were numerically integrated with different values for n 
and b, in units c=l,t p = ^p. 

It has been in particular found that: 

(1) For n — 1, the backward null geodesies starting from (r/ sl ,t lsl ) reach r = 
without leaving the matter dominated area (approximately delimited near r = by 
T = T cq = 10 5 °K) provided b is kept larger than about 10 12 s. 
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(2) for n = 2, an analogous condition holds, the limiting value for b being about 
10 14 s. 

Figures 1 and 2 show the case n = 2, b = 5xl0 14 s for which t c = 2.18 10 12 s. 

5 Discussion and conclusion 

Using a delayed BB universe, the horizon problem has been solved for a class of 
simple models fulfilling some restricting conditions. Further work will be neces- 
sary to discriminate between these conditions which are generic and which are only 
generated by the assumptions made for simplification purpose. 

For instance, the dust approximation used in this letter has been retained to 
allow analytical calculations. It was shown in section 4 that this dust choice provided 
constraints upon the BB function t (r). 

The behaviour of light cones in the radiation dominated region is thus an ap- 
pealing issue for future work. If it could be proven that the geometry of this region 
does not bend the light cones such as to have them reach the BB surface before 
r = 0, then the above constraints could be discarded. 

One could therefore consider the a priori interesting case of a BB function, ful- 
filling conditions summarized at the beginning of section 4, but arbitrarily close to 
the FRW f„(r) = 0. 

In this case, the "unnatural" prescription "centered Earth" would no more be 
needed as the conditions for standard nucleosynthesis could be verified for an ob- 
server located at arbitrary values of r. 
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This would imply that the characteristic lenght of the density inhomogeneities, 
written for the radiation dominated model, at the Earth location, is much smaller 
than the mean free path of the nucleons. A constraint would thus appear on to(r), 
i.e. limiting b to small values compatible with the "close to FRW" assumption. 

The authors are well aware of a potential difficulty of the present model, namely 
to put the observer near the "center" of the universe. In addition to the fact that 
such a location is not forbidden by scientific but only by philosophical principles 
(which they do not accept), they want to stress that the present model is only a 
first "toy-model" . They hope to build, in the future, less simple models, getting rid 
of this prescription. 

In a recent work, they have shown that a delayed BB of type to( r ) — ^ r can 
reproduce the observed dipole and quadrupole in the CMBR anisotropies - a first 
version of this work has been submitted for publication H] with a decreasing BB 
function (negative values for b), but it is easy to see that the same results hold for 
b > 0. 

It comes out from this work that to any given value of the location r p of the 
observer corresponds a value of b for which the observed data are reproduced, b is 
all the smaller as r p is larger. 

A reliable model of universe of this kind could thus get rid of the "centered 
Earth" assumption, provided b should be sufficiently small. This implies that the 
null geodesies, if causally connected, should be so in the radiation dominated region. 

Now, why should we feel unconfortable with the idea that we could be located 



near the "center" of the universe? Following Ellis, Maartens and Nel [|1(| 11], who 
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also dared assume such an "unnatural" prescription in their Static Spherically Sym- 
metric (SSS) model of universe, one can claim that this is no more (un) reasonable 
than the belief in a Cosmological Principle. The purpose is not to put the observer 
at the "center" a priori, but to answer the question as it was put forward by these 
authors in their cited papers: "Given a universe model of the type proposed, where 
would one be likely to find life like that we know on Earth?" . The answer of Ellis 
is: "where conditions are favorable for life of this kind ... near the center, where the 
universe is cool." 

The flat universe approximation, even if it seems more physically justified, can 
also be discussed. If our universe would be proved not so flat as it seems to be, the 
study of the open (closed ?) case would become necessary. 

However, even if the here proposed class of models appears as a restrictive answer 
to the horizon problem, it might, with some easy to conceive improvements, equally 
account for structure formation and all scales anisotropies of the CMBR, which 
could otherwise proceed from a topological defects like mechanism. This will be the 
purpose of other work to come. 

Inflation was, from the beginning of its success story, equally aimed at solving 
the flatness and monopole problems. 

As it is a mere product of Friedmann's equations, the flatness problem only per- 
tains to FRW universes and is thus irrelevant for the class of models here proposed. 

As for the monopole problem, delayed BB without inflation only implies that 
topological defects theories leading to a production of local stable monopoles are 
ruled out. 
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This letter is a first attempt to show that the delayed BB scenario, as it is a 
natural and simple way to solve the problems of standard cosmology while keeping 
its best successfull predictions, is worth spending time and efforts to bring it to at 
least as worthy a paradigm as any other on the market place. 

The last point to emphasize is the following. In the years to come, two satellite 
boarded missions, MAP to be launched by NASA and Planck to be launched by 
ESA, will be dedicated to a high-resolution maping of the CMBR anisotropies. One 
of their main purposes is to provide a test of cosmological theories and an estimation 
of cosmological parameters. 

Number of recent papers attempt to show how the values of these parameters 



could be determined by an analysis of the data thus obtained, see e.g. [fl2[] , |L3 



and []T4 1 . In these papers, the BB function to( r ) is always implicity or explicity set 
to a constant value over the spatial coordinate r, and the cosmological parameters 
considered are those pertaining to a universe with FRW background. 

To be complete, the analysis of these future data will also have to be performed 
in the light of the present results. 
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Figurel. case n = 2, b = 5xl0 14 . The solid lines are the light cones, the upper 
one being the null geodesic issued from (r p ,t p ) ending at (nsij^lsi) on the last scattering 
surface, the lower one being the null cone from (nsij^lsi) to (r = 0,t c ). The dashed line 
represents the Big-Bang surface, and the plotted one, the last scattering and the shell- 
crossing surfaces which cannot be resolved at the scale of the figure. 
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Figure2. Figure 1 zoomed on small values of r and t. The solid line is the light cone 
from (nsij^lsi) to (r = 0,t c ). The dashed line represents the Big-Bang surface, the dotted 
one, the last scattering surface and the dash-dot-dashed one, the shell-crossing and the 
radiation-matter equality surfaces which cannot be resolved at the scale of the figure. 
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